The pion structure function and jet production in 

7P -^ nX 

Michael Klasen |] 

-^ II. Institut fiir Theoretische Physik, Uiiivcrsitat Hamburg, Lurupcr Chaussee 149, 

f^ , D-22761 Hamburg, Germany 

'^ ' Abstract. 

Despite its theoretical and practical importance, the pion structure is still badly 
p<j I constrained, particularly at low Xt^ and in the sea-quark and gluon sectors. Recently 

ZEUS have presented data on dijet photoproduction with a leading neutron, which is 
dominated by slightly off-shell pion exchange and can be used to constrain the pion 
densities down to x^r — 0.01. We compare a recent NLO calculation to the ZEUS 
data and find that the lower gluon densities of SMRS seem to be preferred by the 
^^ , data. Theoretical uncertainties, in particular from the pion flux, are discussed in some 

(<— ^ ' detail. 
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1. Introduction 

As one the simplest QCD bound states and as the Goldstone boson of chiral symmetry 
breaking, the pion is a very intersting theoretical object: Its structure carries important 
implications for the QCD confinement mechanism and the realization of symmetries like 
isospin in nature. It is also of practical importance for the hadronic input to the photon 
structure at low scales. The latter is connected via Vector Meson Dominance to the p 
meson structure, which is poorly known and thus often replaced by the pion structure. 

Unfortunately, determinations of the pion structure have made little progress over 
the last decade. They are based on old Drell-Yan and prompt photon data at fixed 
target energies and large values of the partonic momentum fraction x.,^. Many details 
are still based on pure theoretical assumptions. In order to improve the situation, it 
has been proposed to measure the (virtual) pion structure at low a;^ in deep inelastic 
scattering (DIS) and photoproduction with leading neutrons at HERA [|l|. Since the 
pion is by far the lightest hadron, its exchange dominates the p ^ n transition and it 
will almost be on its mass shell, particularly at small values of the squared momentum 
transfer t' between the proton and the neutron. 

Recently the HI and ZEUS collaborations at HERA have installed forward proton 
and neutron calorimeters and presented first results on the semi-inclusive DIS structure 
functions F2 {/3,Xp,Q'^) and F2 {(3,Xn,Q'^) 0, 01 and on dijet photoproduction 
0, Q. The latter has recently been compared to a next-to-leading order (NLO) 
calculation 0. In Sec. of this paper, we define the dijet photoproduction cross section, 
the kinematic variables, and the experimental conditions of the scattering process. In 
Sec. 1^ we discuss the pion flux and in Sec. ^ the pion structure. The main results are 
contained in Sec. ^ where we compare our NLO dijet calculation with the ZEUS data. 
A brief summary is given in Sec. ^ 

2. Kinematics 

The production of two jets with a leading neutron at HERA is shown schematically in 
Fig. |I]. The corresponding cross section 

depends on the partonic cross section daab^u/dt and on the parton densities in the 
electron and proton 

fa/eiXa, M^) = f' ^/,/,(x„ M^)f,/,{y, Q^), (2) 

Jxa y 



1 d (l-x „) 

xt i Xr. 



h,,{x,,Ml)= / ^ '^h/^{x^,Ml)U/,{l~Xn,f)- (3) 



Both are convolutions of flux factors and parton densities. The photon flux fj/e{y, Q"^) 
is given by the usual Weizsacker- Williams approximation . The ranges of the photon 
momentum fraction y = {P ■ q)/{P ■ k) G [0.2; 0.8] and virtuality Q^ = —q^ < 4GeV^ 



The pion structure function and jet production in 7p -^ nX 

e ► 1 ► e 



Photon Remnant 




Meson Remnant 



P ► 



Figure 1. Schematic diagram for dijet production with a leading neutron in 

photoproduction at HERA. 



are dictated by the ZEUS experiment. The momentum fraction of the leading neutron 

, t') are restricted 

{l-Xn)(m'^-x„m'^) 



j4 



Xn and the momentum transfer t' in the pion flux factor f-jt/pi 

to Xn = {P' ■ k)/{P ■ k) > 400GeV/820GeV and t' = fipr) = 

with pt < Xn ■ 0.66 GeV. The hard jets are defined with the /ct cluster algorithm 

and a maximal distance A. = 



iVi ~ Vj^ + {4>i ~ 'Pj)^ < 1 II]- They are required 
to have transverse energies above Ej^ > 6 GeV and lie in the central rapidity range 
-2 < V^* < 2. 



3. Pion Flux 

In principle, a variety of soft reactions can lead to a leading neutron final state at 
HERA. Four different possibilities are shown in Fig. |^: a) charged pion exchange, b) 
intermediate A production, c) Reggeon (p, 02) exchange, and d) Pomeron exchange. 
Since the momentum transfer is restricted to low values t' < 0.54 GeV^ in the ZEUS 
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Figure 2. Photoproduction with a leading neutron: a) Charged pion exchange, b) 
intermediate A production, c) Reggeon exchange, and d) Pomeron exchange. 
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Figure 3. Photoproduction with a leading neutron in the meson cloud model at 
large (left) and small (right) impact parameter b. 



experiment, the pion as the hghtest particle in the t'-channel dominates. Diagrams, 
where an additional pion is produced in the proton fragmentation region (diagrams b) 
and d)), are experimentally excluded. Reggeon exchange amounts to a small (< 10%) 
cross section increase. 

In the meson cloud model, the photon-proton collision can be viewed as a process, 
where a pionic cloud is stripped off a neutron core (see Fig. ^. As the impact parameter 
b becomes small, additional effects like photon absorption or neutron-rescattering 
may become important. The semi- inclusive ZEUS data show that absorption effects 
decrease the cross section by at most 10 % in the Xn range considered, which cancels 
approximately the increase of the cross section by Reggeon exchange. 

The light cone representation of the proton wave function is given by 



\P) 



\po) + j da;d2pT {<l)^/p{x,pT)J-\p,7r'^) + ^^/p{x,pT)J-\ 



ra,7r+) + 



(4) 



where |po) is the bare proton wave function, given by a three-quark state, and 
(pM/p{x, pt) is the probability amplitude of finding a baryon B (=p,n) with longitudinal 
momentum fraction x and transverse momentum px and a meson M (=7r°, vr"*") in 
the proton. Z is the wave function renormalization constant and is fixed by the 
normalization condition {p\p) = 1. The splitting function or flux factor takes the general 
form 



U/p{l-Xr„t') = J d'^PT\(t>^/p{Xn,PT)\'^S{t' - fipr)) 



(5) 



3Cn9.. 



mrp 



-t' 



It can be calculated in covariant or time-ordered perturbation theory. C„ is the squared 
Clebsch- Gordon isospin coefficient (2/3 for leading neutrons, 1/3 for leading protons), 
l/(m^ — t')^ is the squared pion propagator, and (1 — Xn)^^^"'"'-* ""''"■' accounts for the 
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Table 1. Parameterizations of the parton densities in the pion. They are constrained 
by Drell-Yan (DY) and prompt photon production (7X) data and by the Momentum 
Sum Rule (MSR). 



virtuality and possible reggeization of the pion. The interaction term in the pion-nucleon 
Lagrangean leads to the numerator — t', and off- mass shell effects of the higher |p, vr*^) 
and \n, 7r+) Fock states are modeled by a form factor 

exp[6(t — m^)] [Exponential] 

^ exp[i?2(t -ml)/{l-Xn)\ [Light Cone] 

The flux factor is constrained by a symmetry relation (/7r/p(l — Xn,t') = fn/p{xn,t')) 
and charge and momentum sum rules. Since the momentum transfer t' is small, we 
neglect reggeization (a' = 0) and choose a light cone form factor with R = 0.5 GeV~^ in 



F{Xn,t) 



(7) 



good agreement with recent determinations |^, |TOl . The pion-nucleon coupling constant 
(7^/(47r) = 14.17 is taken from a recent extraction from the Goldberger-Miyazawa-Oehme 



sum rule |TI| 



4. Pion Structure 



Traditionally, the valence quark densities in the pion have been determined in the Drell- 
Yan process, while the gluon density has been constrained in prompt photon production 
(see Tab. 1). In both cases the relevant data sets were taken in fixed target collisions 
in the 1980s and are restricted to large x^r > 0.2. Furthermore, they suffer from nuclear 
effects and fragmentation contributions, respectively, and also from scale uncertainties. 
The situation is worst for the sea quarks where theoretical assumptions vary most 
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Figure 4. Momentum fractions carried by sea-quarks and gluons in the pion as 
determined in the SMRS fit. The lowest x^ values (6-15) are obtained with sea quark 
momentum fractions between 10 and 20 %. Hi 



widely. It is known that at Q^ = 4 GeV^ valence quarks carry about half of the pion's 
momentum. At the same scale gluons carry about half of a proton's momentum. If 
pions and protons had a similar structure, as assumed by GRSc |T^, no room would 
be left for the sea quarks. Indeed GRV assumed a zero-input density for sea quarks, 
however at a very low scale |T^. SMRS determined that sea quarks should however 
carry between 10 and 20 % of the pion's momentum (see Fig. Q) fl^. In Fig. |^ we 
compare the available parton densities in the pion at a scale Q^ = E^ ^^^ = 36 GeV^ 
relevant for dijet photoproduction on a slightly virtual pion. As one can see, the SMRS 
sea quarks and gluons turn out lower than the very similar ABFKW, GRV, and GRSc 



distributions [|T5|] . 



5. Dijet Production 

The ZEUS collaboration have compared their dijet results with a leading neutron with 
LO Monte Carlo predictions and found reasonable agreement. It is, however, well- 
known that pure LO results are strongly scale dependent, that they can be modified by 
potentially large NLO corrections, and that they do not allow for the implementation 
of a jet algorithm. For inclusive dijet photoproduction, four different NLO calculations 
have been performed using different phase space slicing methods and the subtraction 



method. They all agree well within numerical errors |16]. We have implemented in our 
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Figure 5. Parameterizations of the valence-quark (top), sea-quark (center), and 
gluon (bottom) densities in the pion in LO (left) and NLO (right) at Q^ = 36 GcV^. 
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Figure 6. Dependence of the dijet photoproduction cross section with a leading 
neutron on the transverse energy of one of the two jets. Two NLO predictions 
with different pion structure functions and one LO prediction are shown together 
with experimental data from ZEUS with statistical (inner) and systematic (outer) 
experimental error bars. In addition, error bands from the theoretical scale uncertainty 
(medium: LO, dark: NLO) and the experimental energy scale uncertainty (light) are 
shown. The NLO error band coincides with the thickness of the full line. 



calculation [jT^ two additional integrations, i.e. over x„ and t', and the pion flux factor. 
The proton parton densities have been replaced with those in the pion. 

Selecting the two jets with highest Et and integrating over the t) G [—2; 2], we 
compare in Fig. ^ the LO and NLO Et distribution with GS96 photon densities [|18] and 
two different pion densities to the ZEUS data. Almost identical results are obtained with 
GRV photon densities. |19[ The i^-factor turns out to be one, but the scale dependence 
is reduced considerably from LO to NLO. The normalization is sensitive to the chosen 
pion flux factor: If one includes the Regge trajectory [a' = 1 GeV^^) and omits the 
exponential, the cross section is reduced by 15%. The same reduction is obtained with 
a light cone form factor and R = 0.6 GeV~^. Resolved photons and gluons in the pion 
are important at low Et < 9 GeV and Et < 22 GeV, respectively. In the ratio of 
leading neutron and inclusive dijet cross sections in Fig. |^, many uncertainties from the 
scale choice, photon densities, and hadronic energy scale cancel. The Et shape is then 
sensitive to the pion densities and shows a preference for SMRS3 over GRSc. 

If the leading neutron cross section is integrated over the two highest jet Et > 
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Figure 7. Ratio of the leading neutron over the inclusive cross section versus the 
transverse energy. NLO calculations with two different pion structure functions are 
shown together with the experimental ZEUS data. Only statistical error bars are given. 



6 GeV, but not over 77, one obtains the rapidity distribution shown in Fig. ^ Direct 
photons and quarks in the photon are important at 77 < 0. In this region of jet 
photoproduction, the low quarks in the NLO GS96 parameterization usually give the 
best fit, despite the fact that they do not reproduce the photon structure function. [^ 
Indeed GS96 also gives the best fit for dijets with a leading neutron. At 77 < 
hadronization corrections become also important and can amount to a factor of two. 
The data for the ratio of leading neutron over inclusive rapidity distributions in Fig. 
^ have still rather large statistical error bars and do not distinguish between different 
pion parameterizations. 

The distribution in the observed momentum fraction of the partons in the pion 



X 



2fl 



X. 



,E 



(8) 



is shown in Fig. |T0|. It suffers from the fact that ZEUS have applied equal Et cuts on 
both hard jets, which is known to be infrared sensitive [^. Therefore the cut on the 
second jet in the NLO calculation has to be relaxed to Et,2 > 5 or 5.5 GeV. The best 
fit is again obtained with GS96 photon densities and SMRS3 pion densities, which have 
the lowest quark and gluon distributions, respectively. It should be kept in mind that 
the quark densities in the photon are correlated with the gluon density in the pion. In 
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Figure 8. Dependence of the dijet photoproduction cross section with a leading 
neutron on the rapidity of one of the two jets. Details are as in Fig. g. The lower full 
curve demonstrates the influence of hadronization corrections. 



the XTT-distribution, hadronization corrections may play an important role at small Xtt, 
which is related to negative rapidities. 



6. Summary 

Despite their theoretical and practical importance, the parton densities in the pion 
are still poorly constrained, particularly below x.,^ < 0.2. An intriguing possibility for 
improvement consists in dijet photoproduction with a leading neutron at HERA, which 
is dominated by slightly off-shell pion exchange. Here the pion can be probed down 
to Xtt — 0.01. However, theoretical uncertainties from Reggeon exchange and photon 
absorption, modeled by a pion flux factor, have to be well under control. 

Data on dijet photoproduction with a leading neutron have been obtained by 
the ZEUS collaboration. Unfortunately, ZEUS have applied theoretically unsafe Et 
cuts. Furthermore, the data suffer from large hadronization corrections in the backward 
region. On the other hand, the hard scattering cross section has been calculated in NLO 
and is theoretically stable. A comparison of the NLO calculation with data indicates 
that lower gluon densities in the pion, as predicted by the SMRS3 parameterization, are 
preferred over steeper gluons as predicted by ABFKW, GRV, or GRSc. 
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Figure 9. Ratio of the leading neutron over the inclusive cross section versus the 
rapidity. Details are as in Fig. 0. 
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